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Abstract 


A new technique is proposed for measuring volcanic deformations on the 

order of centimeters per day to centimeters per year, called the Airborne Radar 

Volcanic Surveillance System (ARVSS). Such a system would employ an airborne 

multi frequency pulsed radar, tracking passive ground reflectors spaced at 1 

kilometer intervals over a 50 square-kilometer area. Identification of targets 

would be done by doppler and range resolution techniques, with final relative 

position measurements done by phase comparison of multi frequency signals. The 

system is impervious to all but severe weather, and can be operated day or night. 

Atmospheric path length errors may be corrected by an airborne refrac tome ter, 

meteorological instruments, or other refractive measuring devices. Anticipated 

A 

system accuracy is 1-2 cm, with measuring times on the order of minutes. Poten- 
tial problems exist in the high intrinsic data assimilation rate required of 
the system to overcome ground backscatter noise. 

Present geodetic techniques are too slow, cumbersome, and costly to be 
used on many remote volcanoes. Alternative proposed geodetic systems suffer from 
various drawbacks which make them unsuitable as volcanic deformation measuring 
tools. It appears that ARVSS satisfies all the requirements for such a tool. 
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Introduction 

The author worked during summer 1980 as a Visiting Graduate Fellow at Lunar 
and Planetary Institute, on the development of a remote sensing system for detec- 
tion of volcanic surface deformations. As a Summer Intern In 1979, under the gui- 
dance of Drs. Stan Zlsk and Tom McGetchln, the preliminary design of such a system 
was completed. This summer, with Zlsk and Or. Roger Phillips as advisors, an analy- 
sis of data processing requirements for the system, as well as a more definitive 
study of slgnal-to-nolse requirements, were performed. As a result of these analy- 
ses, parameters such as reflector size, signal Integration time, end effective 
bandwidth were determined for various system scenarios, which will form the basis 
for final system selection and hardware construction. Details on the system are 
contained in the draft of a report by Zlsk and Lunine entitled "Airborne Radar 
Volcanic Surveillance System". A digest of the results from the report appears 
below. 

The Problem 

Ground surface deformations associated with volcanic activity are believed 
to be caused by subsurface magma movement. During active periods areas near the 
cone of a volcano may deform at a rate of centimeters per day, with areas on the 
flank moving at a rate of centimeters per year or slower. Measurement of such 
activity, to an accuracy of a centimeter or better, over a 50-100 square-kilometer 
area, is a major concern of a number of volcanologists. Magma chamber shape (Decker 
and Klnoshlta, 1971) and lateral magma transport can be deduced from horizontal 
and vertical deformation data. The possibility also exists that such Information 
can be used to predict eruptive or Intrusive activity. 

Present techniques for measuring volcanic deformations are all ground-based. 
Tlltmeters, while rapid and Inexpensive, cannot give complete information on ver- 
tical and horizontal motion. Optical levelling provides vertical movement data to 
millimeter accuracy, hut requires days to cover a 50 square-kilometer area, and 
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Is too costly to do frequently (Decker and Klnoshlta, 1971). Relative horizontal 
movement Is determined by ground electronic radio-ranging devices, which provide 
better than millimeter accuracy after atmospheric correction. Such ground- level 
devices suffer from llne-of-slght difficulties, are time-consuming for a large 
reflector field, and are sometimes dangerous, as Indicated by the death of a 
geologist performing such measurements on Mt. St. Helens at the time of the May 
18, 1980 eruption. 

Thus only a few volcanoes In Hawaii, Japan, and Iceland have been subjected 
to detailed and frequent deformation studies. A number of geologists have expressed 
the need for a relatively Inexpensive tool for measuring volcanic crustal movements, 
with an accuracy In the range of one centimeter, and data collection times sig- 
nificantly less than a day. A number of current and proposed techniques, both 
ground-based and air/space-borne, have been devised by various groups for more 
general tectonic work, but none fits the requirements for a volcano-oriented 
system, based on analysis by the author. Thus an entirely new system, the Airborne 
Radar Volcanic Surveillance System (ARVSS), was designed by Zlsk and the author 
to satisfy the requirements of volcanologists. 

The System 

ARVSS involves the placement of passive radar retroref lectors on the vol- 
cano of Interest to form a reference grid, from which displacements of the sur- 
face can be detected. To adequately detail the ground movement, the reflectors 
would be spaced at 1 kilometer intervals, and cover roughly a 50 square-kilo- 
meter area (100 square-kilometers in larger rift systems). To achieve a high 
signal return over a broad angular range, the reflectors would be of the trihedral 
corner type, consisting essentially of two perpendicular specular bisecting planes. 
They would be anchored In the ground and covered with a protective foam enclosure 
to prevent coating by volcanic dust and water. 

Remote detection of the passive retroref lectors would be performed by an 
airborne radar device, which would be flown over the volcano a sufficient number 
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of times to detect all of the reflector grid points. Processing of the returned 
radar signal, for range and range rate (doppler) Information, would yield the 
relative reflector positions to within meters, and phase comparison of the multi* 
frequency signals would then provide the fine (1.5 cm) resolution of relative 
reflector positions. Performance of additional flyovers after a time Interval 
ranging from days to months would reveal relative displacements In reflector 
positions, thus providing a direct measure of volcanic ground surface movements. 

The type of radar, reflector size, and form of processing employed to pro- 
vide the centimeter-scale accuracy required have been the subject of study by Zlsk 
and the author. The results are summarized below. A flowchart of the analysis 
performed i$ shown In figure 1. The following are "constants" chosen for the system: 

A system frequency of 10 GHZ (X-band) has been selected for Its short" wavelength 
(hence high Intrinsic accuracy— 1-2 cm), and low water vapor signal attenuation 
(a problem for shorter wavelengths). A medium-sized carrier aircraft Is assumed, based on 
user economic constraint, with a cruising altitude of 3 km. and a velocity of 
0.0894 km. /sec (200 mlles/hr.). The velocity of signal propagation through the 

c 

atmosphere, c. Is set as 3.00 X 10 km. /sec. 

Cell Area Calculation : Consider a square trihedral reflector of panel langth a. 

The cell area on the ground which must be resolved to detect the reflector Is 

C*i<* n/e«r* t !~ ■ — ■ a ' —- (1) 

y s/‘j * 1 

where a«ref lector panel length 

i *wavelength 

S/N s s1gna1-to-no1se threshhold detection level 
r . “ground backscatter per unit area 

Values of A for representative volcanic surfaces. Incidence angles of the radar 
beam, and reflector sizes are shown In table 1. 

Doppler Resolution : ARVSS resolves the ground Into a matrix of cells. Doppler 

resolution fixes cell width In one direction, while antenna beamwldth or, If 
necessary, range resolution constrains the perpendicular direction. A study of 
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what could be achieved with doppler resolution was performed to determine whether 
range resolution would be required to limit cell size. 

Consider the radar aircraft to be at the origin of a Cartesian coordinate 
system, the aircraft moving along the y-di recti on with a velocity -v y . Then a 
reflector, placed at an arbitrary position In the coordinate system, may be seen 
to move with a velocity +v y In a frame Initially coincident with the original 
coordinate system, but moving with the aircraft. This frame Is shown In figure 2. 
Then 

cq ^ 1 - 


,* £ 
Vj f* 




t > ^ ' (2) 


As the aircraft moves along Its path, the projections of the surfaces of constant 
doppler frequency on the ground form, for finite y, a family of proper hyperbolae. 
This can be seen by setting v f , the radial velocity, and z constant. In an 
equation similar to (2) but expressed In terms of v f rather than f^. Me get 

y >( v/- 1 //.) - i ,'// + v/i 3 (3) 


which Is the equation of a proper hyperbola for finite y (Korn & Korn, 1968). 
Hyperbolae of constant doppler are plotted In figure 3. 

The change In doppler frequency with time Is 

c r t ( 4 ) 


A,* 




C (jL^ % X **tV 7 '* 


which Is plotted for various x and y In figure 4 (y Q in the figure Is the Initial 
y-posltlon of the aircraft— y Is continually recalculated as t Increases). 

The doppler resolution cell width q (see figure 5), Is a function of f d and 
the doppler frequency spacIngAf^. Explicitly: 



i 

Jr V s « ex V* 



(5) 


Note that qecAf^. Values of fj and q for Jt f^ B 2 Hz are given In table 2. 

Given the doppler cell widths shown In table 2, Is range resolution necessary 
to achieve a cell size sufficiently small to permit target reflector detection? 

We consider a cell constrained In the range direction by the major axis of the 


projected beam on the ground (a "worst case"). Then 

CGll W31k iri i /f c ^ ^ 

where 4* -antenna beamwldth, $ “beam's angle of Incidence from vertical. Then 

A*PP*q . and from (1): « 

o . - a c o- . 

,/V, jM ^ 

Values of S/N, for given & , <E and q from tables 1 and 2, are given In table 3. 

For a "reasonably-sized" reflector (a£0.5 m), the S/N Is much lower than 
the 15 dB required to detect the reflector. Thus range resolution (via multifre- 
quency or pulsed radar) Is needed to limit the range direction cell width to a 
size much smaller than Is attainable by beamwldth limitation only. A 3. meter 
reflector can be detected within the beamwldth-llmlted cell, and several might 
serve as orienting points for the field, but such a complication will not be con- 
sidered here. 

Range Resolution : To detect the *$neter panel length reflectors requires, from 

equation (1), A=25 m . For q*l m, range cell width s must be 25 m. This requires a signal 
bandwidth f 

8 = sj * ^ ( 8 ) 

Assume a multi frequency CW (continuous wave) system. To avoid ambiguities In the 
signal return, the minimum frequency spacing within B must be 

(9) 

where d>al1as1ng (range ambiguity) distance, chosen to be the length of the reflec- 
tor field (*> 7.5 km). 

Then to achieve the effect of a random signal with the above minimum fre- 
quency spacing (300 such spaclngs spanning B), which Is necessary to detect the 
reflector peak In the returned signal, 24 frequencies are required.* To achieve 
the fine resolution of 1.5 cm, 50 frequencies are required, for each of the 24. This 
Is a prohibitively large number of frequencies from a practical standpoint. The 


*A "random signal" Is created by selecting frequencies within the bandwidth with all 
possible spaclngs, from 20 KHZ up to the bandwidth. This requires 24 frequencies. 
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required number of frequencies can be drastically reduced by pulsing the radar signal. 
Pulsed System Scenarios : Briefly, the processing scheme Is as follows: compressed 

pulse signals, emitted from the airborne transmitter, bounce off the ground and 
are received by the aircraft. These signals are processed to provide rough reso- 
lution of the reflectors from doppler and range Information. This Is done at a 
number of frequencies simultaneously, which are then phase-compared to provide 
fine resolution of the reflectors' positions. 

Three such systems are detailed In table 4. For the purpose of the analysis, 
a projected cell width (for the H meter panel length reflector) of 15 m, and 
projected doppler resolution of 4,2 Khz are considered. The equations used In 

the analysis are outlined below. 

h - 4 ™ 5a***vj. (1Q) 


of - (11) 

if at " f<v!« (^) 

FFT processing equations— To process the Incoming signals, a Fast Fourier Transform 
(FFT) procedure Is used. The number of samples per range box which must be taken 
to unambiguously detect the reflector signal Is 


. (13) 

The FFT Is not straightforward— as the aircraft moves along Its path, reflectors 


ox f *>pi?tssfc ~ 


will shift out of cells (slip) and processing must be done In frequency sections 
to avoid losing the reflector signal. The maximum doppler frequency which can 
be processed without slippage Is 

f* - m* m 

J S A 

and processing must be done In frequency bins of width f ds* from ' f d max 10 +f d ma> 
Finally: jaf. p, c r.v* 

v /csJ. 

Discussion : In each of the cases shown In table 4, recording rate demands are 
excessive, particularly for the 0.1 microsecond pulse system. Current tape devices 
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record at 60,000 megabytes/ second, while fast computer memory storage and optical 
recording devices are currently prohibitively expensive for this application. Beam- 
ing the signals from the aircraft to ground recorders Involves transmission compli- 
cations which make this alternative Impractical. The recording problem has not 
yet been solved. 

The number of frequencies required for fine resolution Is another system 
constraint. Over 3 times as many frequencies are required In the 0.1 microsecond 
pulse case as In the 0.01 microsecond pulse case. Onboard FFT processing to Identify 
cells containing the reflectors ( and hence reduce the number of wells which must 
be processed for fine resolution) Is most easily done In the 3 meter reflector case. 
However, such a large reflector Is not really a practical alternative. For the 
H meter panel length reflector case, the number of samples required for the FFT is 
excessive for the 0.1 microsecond pulse, and somewhat less so for the 0.01 micro- 
second pulse, primarily because of the reduction in the number of required frequen- 
cies. Although recording speed Is very high In the 0.01 microsecond pulse mode, 
real-time processing need not be— one can record the returned pulse in 1/10 of a 
second, and then no pulses need be recorded over the remaining 9/10 of a second. Then 
the same amount of time Is available for processing as In the 0.1 microsecond pulse case 

The problem of scintillation, or noise adding In phase In the reflected 
signal, was also analyzed, and found not to be a problem for the processing 
considered here. 

Atmospheric refraction as a source of reflector position error has been 
examined. There are two components to this error: a stable "dry atmosphere" effect, 
and a highly variable water vapor Influence. The latterftay cause an apparent 
path length extension of up to 25 cm. A number of means of measuring the Index 
of refraction or water vapor content of the air have been reviewed by the author. 

It appears that standard airborne meteorology equipment can correct for the effects 
to within 1-2 cm, and a new airborne refractometer developed In the Soviet Union 
by Levin, et al (1977) may provide correction to a centimeter. 
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Conclusion 

Despite the fact that no clear-cut choice between the above pulse systems 
can as yet be made, and that serious recording rate problems exist In all cases 
considered, the system parameters which have been calculated Indicate that the 
ARVSS type of deformation measurement can be performed. The measurement can be 
„ done day or night. In most weather, within the space of a half-hour. The use of 
. passive reflectors eliminates the need for servicing visits to the observed 
volcano. Additional applications of the system Include mapping of subsidence zones, 
saltdome deformation Investigation for selection of nuclear waste disposal slttf, 
and mapping of extraterrestrial volcanic deformations. The system's primary 
Intended use, however, Is In the measurement of terrestrial volcanic deformations, 
in areas Inaccessable to current measuring devices. 
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Figure f _ 

Flowchart of AR7SS Analysis 


Cell size as a function of 
reflector panel length, 
surface type, vertex angle, 
signal -to-noise ratio thresh 
hold. 


Doppler frequency geometry: 

1. Derivation of radial velocity v f . 

2. Derivation of locus of constant v„. 

r 

3. Plot of doppler frequency as a function 
of time. 

4. Calculation of doppler cell width q 
as a function of doppler frequency, 
doppler frequency spacing, aircraft 

height above ground, azimuth and vertex J 
angles. ! 

5. Cell area given range cell length j 

constrained by antenna beamwldth only. 

Signal -to-noise ratio of resulting cell area, i 
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S1gnal-to-no1se ratio for cell area 
constrained in range direction by 
multi frequency range resolution; number j 
of frequencies required. 

— j 


Construction of possible ARVSS 
configurations. 
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